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Abstract 

, The Brownian dynamics of the density operator for a quantum system interacting 
with a classical heat bath is described using a stochastic, non-linear Liouville equa- 

fSJ ' tion obtained from a variational principle. The rate of the noise-induced transitions 

^ . is expressed as a function of the environmental spectral density, and is discussed 

I for the case of white noise and blackbody radiation. The time-scales of decoherence 

' and dissipation are investigated numerically for a system of two quantum states, 

i These are the ground and first excited state of the center of mass vibrations for an 

, ion confined in a harmonic trap. 
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, I. Introduction 

^ , . The evolution of a quantum system surrounded by a thermal environ- 

ment represents a fundamental problem of the quantum theory, increasingly 
important for applications. In classical mechanics, a particle interacting 
^ ■ with a heat bath has a Brownian evolution towards thermalization, which 

reduces to a purely dissipative motion when the temperature decreases to 
zero pp. For a quantum particle, the interaction with the environment is 
more complex and may produce phenomena without a classical analog, as 
the collapse of the wave function and decoherence P]. These effects are 
particularly important for understanding the mechanism of atom switch- 
ing in the scanning tunneling microscope 0, for the practical realization 
of a quantum computer jlj 0! or the transmission of quantum information 
j^. However, despite the growing interest in these fields, the physics of the 
Brownian quantum dynamics and of the wave function decoherence is not 
yet completely understood [7j |S]. 

The conceptual difficulties appearing in the description of a quantum 
system interacting with a classical environment are well illustrated by the 
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" Schrodinger cat experiment", or the Zeno paradox at the continuous mea- 
surement process jH] ^Oj. The wave function collapse at macroscopic scale 
(e.g. of the apparatus pointer) was investigated by Zurek [2j, assuming that 
quantum mechanics describes both the environment and the system. How- 
ever, when only the system is quantum, the reduction of the density operator 
by partial tracing over the environment's degrees of freedom is not possible, 
and the approach should be different. The problem occurs, for instance, in 
quantum gravity, whenever the space-time framework of a quantum particle 
remains classical 11 . In such cases, a promising approach to the mixed 
classical-quantum dynamics consists of using coupled Hamilton-Heisenberg 
equations derived from a variational principle . The time-dependent vari- 
ational principle (TDVP) jl3j is frequently used in the quantum many-body 
theory to obtain the mean-field equations (e.g. time-dependent Hartree- 
Fock). The mean-field dynamics appears by kinematical constraints on the 
unitary evolution in the many-body Hilbert space, and has a mixed char- 
acter, including both quantum and classical aspects |14j . Therefore, the 
TDVP provides an appropriate frame to treat also the mixed version of 
the bilinear coupling model, consisting of a quantum particle interacting 
with a heat bath of classical harmonic oscillators. For the pure states, this 
approach leads to a non-linear time-dependent Schrodinger equation, con- 
taining additional noise and friction terms [^. However, for the treatment 
of the thermalization process, and evolution from pure to mixed states, it is 
necessary to include the interaction with the environment in the dynamical 
equation for the density matrix |16j . 

In this work, the combined effect produced by noise and dissipation on 
the quantum evolution will be described using a stochastic, non-linear Liou- 
ville equation. This equation can be derived from the variational principle, 
and will be presented in Sec. II. 

The noise represents an external stochastic force which stimulates transi- 
tions between the stationary states of the quantum system. A perturbative 
calculation of the transition rates is presented in Sec. III. The result is 
compared with the one obtained when the environment is quantized, and 
is illustrated by two examples: the Gaussian white- noise and the thermal 
radiation. 

The evolution of the density operator when noise and dissipation are 
both included is investigated numerically in Sec. IV, for the case of a quan- 
tum two-state system with relevance for the operation of a logical gate in 
a quantum computer. This system consists of the ground and first excited 
state of the center of mass vibrations for a cooled ion in a harmonic trap. 

The summary of the results and the conclusions are presented in Sec. V. 
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II. The stochastic Liouville equation 



Let us consider a quantum system with the Hamiltonian and density op- 
erators Hq and p, respectively. These operators are Hermitian, and their ma- 
trices in an orthonormal basis {|e„)} can be diagonalized by unitary trans- 
formations. If denotes a normahzed eigenstate of p{t) corresponding 
to the eigenvalue Wn, then there exists an unitary operator Xi such that^ 
p = x^x^ ^ with D = J2n''^n\en){en\- The operator D is related to the von 
Neumann entropy S = —Tr{plnp) = —Tr{DlnD) = —J2n''^n^^Wn- 
According to Reznik ^Hl it is convenient to introduce a "square root oper- 
ator" T], such that p = rji]^ . This operator can be defined as r] = xVD^ and 
in this case D = ri^rj = p — [rj, r/^]. 

Following the assumptions of the bilinear coupling model jl9| |2Uj . a 
thermal environment can be simulated by a classical heat bath of harmonic 
oscillators. These will be supposed to be continuously distributed in fre- 
quency with the density /i^^, and interacting with the quantum system by 
the Hamiltonian Hcoup = K dujpujCujQuj- Here K denotes a Hermitian 
coupling operator, are constants, and are the time-dependent bath 
coordinates. The evolution of this mixed classical-quantum system will be 
described using the variational equation 

//•oo 
di{j^ dujpu^{qu,Pu, - h^) + Tr[rjUhdtr] - r]\Ho + Hcoup)r]]} = , (1) 

where the dynamical variables of the quantum system are r], rj^ , and h^j = 
{pIj + m^uP'q^^l^m^^ is the classical Hamiltonian for a bath oscillator. The 
equations of motion obtained from Eq. (1) at the variation of i]^{t), 
q^{t) and pt^{t) are 

roo poo 

ihdtr] = {Ho + K I du:p^C^q^)ri , ifidtri^ = -r]\Ho + K du:p^C^^q^^) 
Jo Jo 

(2) 

for the quantum system, and 

(juj = — , Piv = -iriuiuj'^qu; - CujTr{Kp) (3) 

for the classical oscillators. The system of Eq. (2) leads to the Liouville 
equation for the density operator, 

indtp=[Ho + K dup^C^q^,p] . (4) 
Jo 



^In the finite-dimensional case we may consider also the Gauss decomposition p — 
xDx^ where the matrices % form a nilpotent group, and D is diagonal. 
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The classical equations (3) can be solved in terms of the unknown function 
of time Qp{t) = Tr{Kp{t)), and their retarded solution is 

Qui{t) = [qui{0) + gu>Qp{0)] cosiot + Si^^i^ sin Ljt 

to 

~ QujQ p{t) + Quj I dt'Qp{t') COS iv {t -t') , (5) 
Jo 

with g^) = Cuj/{m^oj'^). When this solution is inserted in Eq. (4), the 
Liouville equation becomes 

ihdtp=[Ho + Wit),p] (6) 

where W{t) = -K[C{t) + /^(t) + Qp(t)r(0) - Qp(0)r(t)]. The additional 
term W contains the force functions ^{t) and fp{t), 

= - r dup^C^ [q^ (0) cos ut + sin ut] , (7) 

Jo 

/■* . 2 r°° ens ii}f 
fp{t) = - V{t-t')Qp{t')dt' , m = - dwJ{u) , (8) 

Jo -K Jo OJ 

where J{u)) = T:uigi^C^p^/2 is the spectral density of the environment. The 
operator K[Qp{t)V{{)) — Qp{'d)V{t)\ ensures the invariance of W to the renor- 
malization of F by an additive constant, and vanishes when r(t) ~ 5{t) or 
Qp{t) = 0. In the case of a classical heat bath at thermal equilibrium we 
will assume that Qp and F are such that Qp(t)T{0) — Qp{0)T(t) = 0. 

The force function fp makes Eq. (6) non-linear, and therefore, in gen- 
eral, the solutions p{t) will not satisfy the superposition principle. Though, 
exceptions may occur on distinguished "superselection sectors", as the triv- 
ial one defined by the set {p/Qp(t) = 0, i G [0, oo)}. 

If the initial coordinates and momenta qu,{0), P(^{0) are distributed within 
a statistical ensemble £ with the temperature T, then ^(t) behaves like a 
stochastic force which has zero mean and is related to the memory function 
r(t) by the fluctuation-dissipation theorem (FDT). Denoting by <C * » the 
average over S, these properties of ^ are expressed by ^ ^{t) S>= and 
^ ^= kBTT{t — t'), respectively. In this case, Eq. (6) becomes a 

stochastic non-linear Liouville equation for the density operator, 

indtp=[Ho-K{m + fpit)),p] . (9) 

Along the Brownian trajectories p{t) obtained by solving this equation, 
Tr{p^) is a constant for every power k = 1,2..., and therefore the den- 
sity operator has a Brownian evolution which preserves the "purity" of the 
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initial state. However, decoherence may appear for the average density op- 
erator pav{t) =^ p ^ {t) = J2r=i Nt P^{^)/^t calculated over an ensemble of 
Nt trajectories p'^{t) generated with the same initial condition, p'^\t=Q = Po- 
This average (supposed to be equivalent to the average over £) accounts for 
the dynamics of the occupation probabilities, and should satisfy a quantum 
transport equation pi OHlinil^ . 

III. The noise-induced transitions 

When the non-linear friction force fp vanishes or it can be neglected, Eq. 
(9) becomes 

ihdtp = [Ho + Hnoise(.t),p] , (10) 

where Hnoise{i) = ~Ki{t) is the noise term due to the fluctuating environ- 
mental forces. This equation has the solution 

p{t) = e-'"^'/^p{t)e'^°''^ , =Te^^o^*'«(*')^i?(*') (11) 
where T denotes the time-ordering operator j23j . 

K{t) = e'^ot/hj^^-iHot/h ^^2) 

is the coupling operator in the interaction representation, and Ca is the 
Lie derivative with respect to the operator A defined by the commutator, 
CaB = [A, B]. The ensemble average of p{t) can be written as 

/5a.(t) = e-*^«*/^/5,,(t)e'^«*/^ , (13) 

and the terms appearing in the expansion of the time-ordered exponential 
defining Pav{t) can be calculated using the FDT. Retaining only the first 
non-vanishing average, the result is 

Pav{t) = PO - ^ dh dt2T{ti-t2)[K{h),[K{t2),Po]] . (14) 

n Jo Jo 

This formula can be used to estimate the rate of the noise- induced transitions 
between the energy eigenstates HE'^)} defined by HqIE/^) = Eif.\Ek). If 
initially the system is in the pure state \Ei), then po = \Ei){Ei\, and the 
rate of the transition \Ei) \Ef) can be defined by the asymptotic time- 
derivative ^ 

A/j = -—^\t^oo (15) 
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of the occupation probability Vf{t) = {Ef\pav{t)\Ef) = {E f\pav{t)\E f) . Us- 
ing Eq. (14), this probabihty is 

Vf{t) = 2A;bT^^^ f dh r dt2T{ti-t2)cos[Qfi{ti-t2)] (16) 
n Jo Jo 

where Kfi = {Ef\K\Ei) is the matrix element of the coupling operator and 
^fi = {Ef — Ei)/h . The memory function appearing here can be expressed 
in terms of J{uj) using Eq. (8), and the transition rate takes the form 

XT'" = ^iKM'kBT^^ , ^/. = |%|. (17) 

This result was obtained assuming a classical environment, but such 
cases are rare. Therefore, it is interesting to compare Eq. (17) with the av- 
erage transition rate A^^ determined by the whole coupling interaction Hcoup 
when the bath oscillators are quantized. Let us denote by \n) the eigenstates 
of Hamiltonian operator /i^ and by p^j the related density operator. At ther- 
mal equilibrium, p^ = Z[J^ exp{—h^/kBT) with = Tr[eyip{—h^/kBT)]. 
The correspondent of the trajectory average pav is Pavi'^) = Tr£[Tl{t)], de- 
fined by the partial trace over the environmental Hilbert space of the density 
operator TZ for the whole system, and within the same approximation as 
above, 

Vf{t) = -^ t dti r dt2{Ef\Tr£{[H,oup{ti),[Hcoup{t2),no]]}\Ef) . (18) 
n Jo Jo 

The interaction representation for Hcoup is defined here with respect to the 
total Hamiltonian Tio = Hq + Hg, Hg = dup^h^^, and 

TZo = \Ei){Ei\(g) Z^^expi-Hs/ksT) , (19) 

with Z£ the partition function of the environment. For an excitation process 
\Ei) \Ef), Eq. (15) with vj of Eq. (18) gives the rate 



2tt °° 

X^Q = -^\Kf,\^ dL0p^y2{p^)nnCl\iqu.)n-l,n\^Sii0fi-u;) , (20) 

^0 n=l 



^2 



consisting of the ensemble average of the Fermi's golden rule [21] • This 
average can be written in the compact form 

X^<^ = Afi{n)^^^ , Afi = ^\Kf,\^J{LOf,) , (21) 
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where Afi denotes the rate of the (non-thermal) spontaneous decay, and 
{n')uj = l/(e^'^/'^sT _ jg |-jjg average number of the bath phonons with 
energy hu). Similarly, the total rate for a decay process \Ef) —>■ \Ei) is A^*^ = 
Afi{{n)i^^- + 1) = X^'^ eyip{hujfi/kBT). At temperatures T » hujfi/kB 
the factor (n)^^. becomes ~ ksT /hufi, and the rates A-*^*^, A^*^ practically 
coincide with xjf"'' of Eq. (17). 

According to Eqs. (8) and (17), the friction force and the transition 
rate are completely determined by the spectral density J {to). In the case of 
Ohmic dissipation J{uj) = 70;, and the transition rate of Eq. (17) becomes 

A^, = ^\Kf,\''^kBT . (22) 

The memory function defined by Eq. (8) is singular, reducing to a delta 
function, T[t) = 2^5(1), and therefore the calculation of the friction force 
requires a special discussion. For a classical particle with the velocity Q, the 
singularity is removed assuming that the cutoff frequency ujc of fii^ is such 
that QjijJc is very small, and the delta function is an approximation for a 
smooth function ^(t) 19 , peaked at t = but non- vanishing for \t\ ^ 0. 
The friction force obtained from Eq. (8) using these assumptions is 

fp{t) = -2^Qp{t) f at - t')dt' = —iQp{t) . (23) 





If 5{t) is supposed to represent the limit of a smooth function defined over 
the interval (— oo,0], then the friction force increases by a factor of 2 with 
respect to the estimate of Eq. (23), being too strong to ensure thermalization 
in a classical system. The power dissipated by the friction force of Eq. (23) 
at T = 0, when the noise disappears, is^ 

^ - -^^^^ = 7«3,f = ^\Trip{H,.m' . (24) 

Though convenient in calculations, the linear dependence of J{oo) on 
uj over the whole spectrum is unrealistic, and it is useful to study also a 
physical situation of interest. Among the various bilinear coupling systems 
[25| ■ particularly important is the case of a charged particle interacting with 
the electromagnetic field. Decoherence produced by spontaneous emission of 
photons from cold two-level atoms during tunneling through a laser-induced 
potential barrier was analysed in ^25^, while a quasiclassical treatment of 

^For comparison, the power dissipated in the quantized environment, by a spontaneous 
decay \Ef) ^ \E,), is Af^{Ef - Ei) = 2j\{Ef\[Ho, K]\E,}f /h'' . 
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the bremsstrahlung spectrum for a tunneling charge is presented in [27j . 
Let us consider here the case of a quantum charged particle interacting with 
the blackbody radiation field. In the dipolar approximation, the interaction 
Hamiltonian has the form H'j^g^ = —d ■ E = —dxE^ — dyEy — dzE^ where 
d = ef is the electric dipole operator of the particle, and E is the radiation 
electric field. For isotropic radiation eE represents a stochastic force with 
mean, such that -ff^of^e ^ sum of three terms, each of the form —K^ 
assumed above. If the components of the stochastic force are uncorrelated, 
the FDT is expressed by < Ei{t)Ej{t') >= 6ijkBTT{t - t'). For t = t', 
the left-hand side of this equality can be related to the total spectral energy 
density u^j, 

<C Ef{t) >= — / duju^ , (25) 

o Jo 

and if T{t) in the FDT has the form assumed in Eq. (8), then 

27r2e2 

For Uuj the realistic expression = tiuj'^ {n) ^ / -k"^ can be considered, which 
accounts for the effects of the second quantization of the radiation field. In 
this case, the transition rate provided by Eq. (17) is 

A;^f = Ay,(n)^,, (27) 

where 

denotes the Einstein coefficient for spontaneous emission. Therefore, A^f IS 
the usual rate of the stimulated transitions in thermal radiation field [28j . 
The memory function vanishes when T = 0, but at high temperatures it 
becomes proportional to the second derivative of the delta function, T{t) = 
— 4e^(5(t)/(3c^), and the corresponding friction force 

is the classical radiation reaction [29j . 

IV. Numerical results for a quantum two-state system 

The model of a quantum two-state system (TSS) interacting with the en- 
vironment appears in many physical situations [SOI- Usually it is associated 
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with the relaxation of a 1/2 spin system, but it was adapted also to study 
the environmental effects on the quantum coherence oscillations (QCO) in 
one-dimensional double-well potentials. 

A quantum TSS can store one bit of information (qubit), and represents 
the main component in the design of a quantum computer. The influence 
of the environment is crucial in this case, because an accurate execution of 
the logical operations pT| requires a conditional, but unitary evolution of 
the qubit states. According to Cirac and Zoller, |32j a first generation of 
quantum logic circuits could be created using cold ions confined in a linear 
trap. The qubits are reprezented by two- level systems of internal electronic 
states, while the conditional dynamics required by the logical operations is 
implemented by their entanglement using the collective states of the center 
of mass (CM) vibrations. Therefore, the operation time of such computing 
devices is limited essentially by the occurrence of spontaneous transitions 
between the electronic states, and the environmental decoherence of the CM 
vibrational dynamics f5H|. 

The present numerical estimates concern the quantum TSS consisting 
of the ground (|0)) and the first excited state (|1)), for the CM vibrations 
of a single ion trapped in a one dimensional harmonic oscillator potential. 
Within this restricted two-state Hilbert space, the unperturbed Hamiltonian 
and the coupling operator can be expressed in terms of the Pauli spin ma- 
trices a, = |1)(1| - |0)(0|, Gy = |1)(0| + |0)(1|, a, = i(|0)(l| - |1)(0|). This 
representation is chosen such that Hq = A(7a;/2, and K = Qaz, simulating 
a coupling linear in the CM coordinate or momentum. For simplicity, the 
spectrum of the environmental noise will be considered flat (white noise), 
having the density J{uj) = 70;. In this case Eq. (9) becomes 

ihdtp=[Ho-QHtyz,p] , (30) 

with J-{t) = £,{t)—^AQPy/h. Denoting by P = {Px,Py, Pz} the polarization 
vector, the density operator can be written as p = {I + P ■ a)/2 and Eq. 
(30) takes the explicit form 

flPx = 2QJ^{t)Py 

nPy = -APz-2QJ^{t)Px (31) 

nPz = APy 

It can be easily seen that P • P is a constant of motion, and along each 
trajectory P(t) the purity of the states is preserved. However, the evolution 
from pure to mixed states may appear for pav =<C p ^= {1+ <C P •f?)/2 
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and is described by the dynamics of the average polarization vector <^ P 
The eigenvalues i^o , wi of ^ p » are related to ^ P » by 



Wo = 



1 + I < P > 
2 



Wl = 



1 - I < P > 
2 



(32) 



and determine the entropy 



S 



{wo In Wo + Wl In wi ) 



(33) 



In the numerical calculations A = hu = 19.74 10~ eV was considered, 
corresponding to a trap frequency w = 30 MHz j34j . The environmental 
decoherence of the CM wave functions in a Cirac-Zoller logical gate is due 
to the heating of the ion vibrational motion, or to the random phase fluctu- 
ations of the laser fields |33] • To simulate the global effect of these external 
factors an effective coupling parameter a = 'jQ'^/h = 10~^ was chosen. 
This value is somehow large, because it corresponds to a spontaneous decay 
rate (Eq. (21)), Aji = 2auj = 6 kHz, greater than the values of practical in- 
terest However, it is small enough to be relevant for understanding the 
behavior of the system at weak damping. Further decrease of a is not ex- 
pected to bring new qualitative features, though the numerical calculations 
may become difficult by the increase of the computer time. The temperature 
was fixed at T = 1 mK, when the thermal energy /c^T = 86.2 10~^ eV is rel- 
atively high compared to A, and the corrections due to the quantization of 
the environment's degrees of freedom may be considered small. The choice 
of a high temperature makes the two-state approximation too restrictive for 
a complete description of the vibrational dynamics, but despite this aspect, 
it remains suggestive because a linear coupling in the CM coordinate or mo- 
mentum produces transitions only between consecutive levels. 

The numerical integration of Eq. (31) was performed using the D02BAF 
routine of the NAG library [SH] using a time step dt = 0.658 ns. The initial 
condition P(t = 0) = Pq was chosen to represent a pure state, and the aver- 
age of the density operator was calculated over an ensemble of Nt = 1000 tra- 
jectories. For each trajectory the noise ^(t) at the moment t„ = ndt was ex- 
pressed by £,{tn) = RnV'^^BT'y /dt where {i?n, n = 1, 2, 3, ....} is a sequence 
of Gaussian pseudo-random numbers with mean and variance 1. This 
choice ensures the FDT in discrete form, <C C(*i)C(*fc) ^= '^ksT^dt^tk/d't. 

The vector Pq was taken of the form Po = (cos 0, sin $) with $ = 
and 7r/2. These values of <I> correspond to an initial preparation of the 
system in the upper stationary state |1), and in the non-stationary linear 
superposition (|1) -|- i|0))/\/2, respectively. 
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Fig. 1. Ensemble average of the polarization vector (A)-(C), and the entropy (D), 
as a function of time. The initial condition Pq = (1,0,0) corresponds to the upper 
stationary state |1). 



The evolution of the average polarization vector and entropy obtained 
when $ = are presented in Fig. 1. The average <^ Px ^ (Fig- 1 (A), solid 
line) is very close to the exponential exp(— 2At) (Fig. 1 (A), dashed line) 
with A = 2akBT/h = 26.2 kHz given by Eq. (22). This result is natural, 
because if the the occupation probabilities vq = {0\pav\0) and vi = {l\pav\i) 
depend on time according to the "phenomenological" rate equations 

Vq = X^vi — X^vq , vi = — A-'-^^i + X^Vq (34) 

then when A^^ A-l- «i A, the relaxation rate of ^ Px f i — vq is 2A. The 
same rate characterizes the decrease of the energy E = Tr{pavHo) = A <C 
Px > /2. 

The averages <C <S Pz ^ presented in Fig. 1 (B),(C), fluctuate 

around zero, and the main contribution to the variation of the entropy (Fig. 
1 (D), solid line) is due to <^ Px The final value of S is very close to 
In 2, (Fig. 1(D), dashed line), the result expected when there is complete 
decoherence and wq = wi = 1/2. Therefore, the pure state |1) evolves to an 
incoherent mixture of the two states, |0) and |1). It is interesting to note that 
the dissipative friction force has a relatively small contribution to the energy 
decrease. This aspect becomes more clear when $ = 7r/2. In this case. 
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t (ms) t (ms) 

Fig. 2. Ensemble average of the polarization vector (A)-(C), and the entropy 
(D), as a function of time. The initial condition Pq = (0,0, 1) corresponds to the 
non-stationary superposition (|1) + i|0))/\/2. 

without environment coupling (7 = 0), Px is a constant, while Py^z perform 
"QCO", Py{t) = sinLot, Pz{t) = cosLot. With coupling, the evolution of 
the average polarization vector and of the entropy are presented in Fig. 2. 
The average ^ Pa; ^ (Fig- 2 (A)) has small fluctuations near the initial 
value, while the oscillation amplitudes of ^ ^ and <^ Pz ^ decrease 
exponentially in time'^ (Fig. 2 (B), (C), dashed line) and asymptotically the 
coherence is lost (Fig. 2 (D)). 

In the calculations presented above the average excitation energy = 
A/2 + Tr{pavHQ) changes in time due to the combined effect of noise and 
friction. When T > a well-defined separation between these two contribu- 
tions is not possible due to the non-linearity, but a measure of the energy 
which is dissipated by the friction force alone is the ensemble average of 
E^iss from Eq. (24). For the TSS considered above this average has the 
form ^ Ediss ^= 0.5A <C P^ where <C Pj' ^ is the ensemble average 

^due to the noise, because for a non-stationary initial state Eq. (14) yields a transport 
equation, in our case with the solution <C ^ <C -Py S>= e~^*e'^*, Q, = V<^^ — A^. 
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Fig. 3. The ensemble average <C » as a function of time when Pq = (1)0,0) 
(A) and Pq = (0,0,1) (B, soUd hne). The short-time behavior of calculated 
without noise (B, dashed line). 



of the component obtained by integrating the equation 



dPl 

dt 



(35) 



The averages ^ P^ ^ corresponding to the trajectories presented in 
Fig. 1 and 2 are pictured in Fig. 3 (A) and (B), respectively, by solid line. 
After a short transitory time, ^ k S> becomes practically a constant Kg ~ 2 
kHz, independently of the initial condition. This constant may provide a 
measure of the effective coupling between the qTiantum system and the ther- 
mal environment, because at equilibrium the average energy dissipated by 
friction should be equal to the average energy which is transferred to the 
system by the thermal noise. 

The short time behavior of <^ P^ ^ is determined by the frictional 
forces and the initial state. If $ = 7r/2 and Eqs. (31) are integrated 
considering ^(t) = 0, then for small times P^ increases linearly in time 
(Fig. 3 (B), dashed line), with a slope kq = = 3 kHz. Asymp- 

totically this linear dependence changes to exponential. The numerical 
calculations show that for t > ts ^ 0.4 ms, P^{t) « Vl{t), Vs{t) = 
n — [n — P^{ts)]eyiY>[—Afi{t — tg)], and the QCO oscillation amplitude is 
aqcoit) ~ aqcoits) exp[-^/i(f - ts)/2]. 

When $ = the friction force vanishes, and without noise {i{t) = 0), the 
system remains in the upper stationary state, without dissipation. However, 
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if $ is very small, but not 0, then agco increases, attains the maximum, and 
then decreases, such that after a transitory time tg, has the exponential 
behavior noticed above, -P^(t) ~ Vl{t). 

V. Summary and Conclusions 

A quantum system in interaction with the surrounding environment has 
a complex dynamics, which cannot be described only by linear unitary trans- 
formations in the Hilbert space between pure states. When the environment 
is classical, the dynamical equations can be obtained using a timc-dcpcndent 
variational principle. In this work the environment's degrees of freedom have 
been simulated by classical harmonic oscillators, while the dynamical vari- 
ables of the quantum system are two non-hermitian " square root operators" 
77, ri\ defined by a Gauss-like decomposition of the density matrix. 

If the coupling is bilinear, the evolution of the density operator is de- 
scribed by a Liouville equation (Eq. (6)) with an effective Hamiltonian 
containing two additional terms. One of these terms (~ ^) is due to the 
external forces, while the other (~ fp) to the retarded backreaction of the 
environment on the quantum system, and is non-linear. Therefore, the back- 
reaction term preserves the purity of the states, but may restrict the validity 
of the superposition principle. 

When the environmental oscillators are distributed within a statistical 
ensemble, the two additional force terms correspond to the noise and fric- 
tion, being related by the fluctuation-dissipation theorem. In this case, the 
density operator has a Brownian trajectory corresponding to the evolution 
from pure to mixed states, described by the ensemble average. 

The ensemble average involved in the evaluation of the transition rate 
between eigenstates can be obtained analytically (Eq. (17)) if the noise is 
treated as a perturbation and the friction is neglected. For comparison, this 
average was calculated also when the environment is described by the quan- 
tum, rather than classical statistical mechanics, using the partial tracing in 
the whole Hilbert space (Eq. (21)). As expected, the two averages become 
close when the thermal energy is greater than the transition energy. 

The perturbative transition rates and the frictional forces are given ex- 
plicitly for an Ohmic environment and for the physical situation of the black- 
body radiation surrounding a quantum charged particle. This latter example 
shows that, by a suitable choice of the environmental spectral density, it is 
possible to recover the stimulated transitions rate provided by the Fermi's 
golden rule. 
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A non-perturbative solution of the stochastic non-hnear Liouville equa- 
tion was obtained numericahy for the case of a two-state system coupled 
to an Ohmic environment. The two states correspond to the ground and 
first excited state of the CM vibrations for an ion confined in a harmonic 
trap, at a temperature of 1 mK. Two different initial conditions have been 
considered, one corresponding to the upper eigenstate of Hq, ($ = 0), and 
the other to an eigenstate of the "coordinate" operator K, = 7r/2). The 
upper energy eigenstate is also a linear superposition between the two "lo- 
calized" eigenstates of K, separated by Ak = 2Q. For this state, complete 
decoherence appears when the excitation energy becomes half of the initial 
value (Fig. 1), after a time ~ td = (2A)~"^ = /{"yksTA'j^), in agreement 
with the previous estimates The decoherence time of the "localized" 
non-stationary superposition is greater by a factor of 2, but the excitation 
energy in this case fluctuates near its initial value (Fig. 2). 

At high temperatures the power dissipated by friction is relatively small 
(Fig. 3), with an equilibrium value which is practically independent of the 
initial state. Without noise, the classical environment produces dissipation 
only if the initial state is non-stationary (e.g. <I> > 0). For asymptotic times 
this frictional dissipation resembles closely the spontaneous decay obtained 
when the environment is quantized. 
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